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Gamma-Ray Bursts: the Isotropic-Equivalent-Energy Function and 
the Cosmic Formation Rate 
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ABSTRACT 

Gamma-ray bursts (GRBs) are brief but intense emission of soft y-rays, mostly lasting from 
a few seconds to a few thousand seconds. For such kind of high energy transients, their 
isotropic-equivalent-energy CEi S0 ) function may be more scientifically meaningful when com- 
pared with GRB isotropic-equivalent-luminosity function (Li so ), as the traditional luminosity 
function refers to steady emission much longer than a few thousand seconds. In this work 
we for the first time construct the isotropic-equivalent-energy function for a sample of 95 
bursts with measured redshifts (z) and find an excess of high-z GRBs. Assuming that the 
excess is caused by a GRB luminosity function evolution in a power-law form, we find a cos- 
mic evolution of E lso cc (1 + z) lm -° a , which is comparable to that between L; so and z, i.e., 
hso 00 (1 + z) 2 - 30 -0 5i (both lcr). The evolution-removed isotropic-equivalent-energy function 
can be reasonably fitted by a broken power-law, in which the dim and bright segments are 
i//(E iso ) oc £-0.27±0.0i and t[f(E iso ) cc E.-0- 87±a07 , respectively (lcr). For the cosmic GRB for- 
mation rate, it increases quickly in the region of < z < 1, and roughly keeps constant for 
1 ^ z < 4, and finally falls with a power index of -3.80 + 2.16 for z > 4, in good agreement 
with the observed cosmic star formation rate so far. 
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1 INTRODUCTION 

Gamma-ray bursts (GRBs) are among the brightest cosmological 
explosions in the universe, mostly lasting from a few seconds to a 
few thousand seconds in soft y-ray. Thanks to quick follow-up ob- 
servations in optical band, the redshifts of some GRBs have been 
measured by detecting the absorbtion lines of their afterglows or 
the emission lines of their host galaxies. So far the number of 
Swi/f-detected GRBs with known redshifts has grown up to about 
one hundred and thus makes a reliable statistical analysis possible. 
Among various statistical works, the luminosity function as well as 
the cosmic formation rate of GRBs are particularly interesting. The 
luminosity function is a measure of the number of bursts per unit 
luminosity, which sheds light on the energy release and emission 
mechanism of GRBs. The cosmic formation rate is a measure of 
the number of events per comoving volume and time, which can 
help us understand the production of GRBs in various stages of the 
universe. 

The isotropic-equivalent luminosity (L iso ) function of GRBs 
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was firstly a ssumed to be a standard candle dWiiersetal1ll998l: 
iTotanil Il999l) and later mo re realistic shapes o f the luminos- 
ity function were derived dFirmani et al. 20041: Guetta & Piranl 



Guetta, Piran & Waxman et alj 120051 : lliang et alj |2007|; 



2005; 

Virgili. Liang. & Zhand2009l) . The cosm ic GRB formation rate has 



also been extensively investigated (e.g., Fenimore & Ramire Ruiz 


|2000j;lNorris et al.l2000j;ISchaefer et al.l200ll: 


Llovd-Ronning et al. 


|2002j; iMurakami et al.ll2003l; lYonetoku et al 


20041; iFirmani et al. 


|2004|; IWanderman & Piranl l20ld: bin et al. 


2010f). For example. 



iFenimore & Ramire Ruia (|2000) found a correlation between the 
variability degree of the prompt gamma-ray light curve and 
the luminosity and then adopted it to estimate the luminosi- 
ties/redshifts of 220 b right l ong GRBs detected by CGRO/BATSE. 
iLlovd -Ronning et alj 12002) used this GRB sample to estimate the 
luminosity function evolution and the cosmic formation rate of 
GRBs. Wei & Gao (2003) find that there is a tight correlation be- 
tween the peak energy of the prompt emissio n spectrum (E p!:iik ) and 
the luminosity. Recently, IZhang et alj d2012h find that both short 
and lo ng GRBs all comply with this correlation. Using this correla- 
tion, |%netota|eral] 12004) estimated the luminosities/redshifts of 
the 689 BATSE GRBs and hence derived their luminosity function 
and the formation rate. Essentially these two works with simulated 
GRB redshifts reached quite similar results, i.e., the GRB formation 
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rate increases quickly in the region of < z < 1 and keeps constant 
up to z ~ 10, which is inconsistent with the cosmic star formation 
rates (SFR s) inferred from UV, optical, and infrared observational 
data so far iMadau et all 1993: iLillv et al.ll99^ : iBarger et al.l2000l : 
IStanwav et alj|2003h . 

The original concept of luminosity function comes from astro- 
physical objects such as stars and galaxies which are long-lasting 
and quite stable in releasing their energy. For GRB-like high en- 
ergy transients, the total isotropic-equivalent energy (£j so ) released 
in the whole duration of one event can be reliably measured, and its 
function (i.e., the number density of bursts per £ iso interval) likely 
provides an independent or even more representative clue on the 
underlying physics. That's why in this work we focus on the so- 
called "isotropic-equivalent-energy function" rather than the tradi- 
tional luminosity function. 

This paper is arranged as follows. § 2 introduces our sam- 
ple and data selection. § 3 presents the statistical technique while 
§ 4 shows the results. We adopt a robust, nonparametric statistical 
technique to derive the isotropic-equivalent-energy function and the 
cosmic formation rate of GRBs from a E lso —z sample. For compari- 
son, the results from a Lj so -z GRB sample are also presented. In §5, 
we discuss the implication of our results and compare the cosmic 
GRB formation rate with the observational cosmic star formation 
rate. Throughout the paper, we use the standard A cold dark matter 
cosmology with the typical parameters fl m = 0.27, Q A = 0.73, and 
h = 0.7. 



2 DATA ANALYSIS 

In this work , two sets of data are ana lyzed. The - z sample 
comes from I Amati et al. I fcOOSt l2009h . containing 95 long GRBs 
and X-Ray Flashes (XRF, i.e. particularly soft bursts). This sample 
is m ade up of two parts . The first part consists of 70 long GRBs 
from lAmati et al.1 J2008h and their redshifts range from 0.033 to 
6. 3. The E iso values in this work are slightly different from those 
in lAmati et alj d2008h because of the different cosmolo gical param- 
eters a dopted in two works. The second part is from I Amati et al.l 
d2009l) without any modification. All the GRB spectra have been 
extrapolated and corrected to [1,10000] keV in the cosmological 

restframe. 

The Li s0 - z sample is from Wanderm an & Pirar] d2010h (Lj S0 
represents the isotropic peak luminosity, and the time resolution is 1 
s). Due to the Swift/BAT narrow energy band, only a small fraction 
of bursts have a well determined spec trum. In order to obtain rea- 
sonable estimates of Lj so for all bursts, Wanderm an & Pirar] j201Ch 
considered the characteristic Band function, i.e., 



N(E) 



M- 



(a-P)Eo 1g - 
lOOkeV 



for£ < (a-p)E 
f^expOS - or), for E > (a - p)E 

(1) 

where the characteristic parameters (E pca ^, a, ji) taken as 
(511 keV, - 1, - 2.25), respectively. To account for the so-called 
k-correction, again all spectra have been extrapolated and corrected 
to [1,10000] keV in the cosmological restfram e. To test whether the 
above Band function applicable to all bursts. Iwanderman & Pirar] 
< |2010|) preformed Monte-Carlo simulation in order to compare 
GRBs having simulated spectraal parameters with GRBs having 
measured spectral parameters. The result of such a simulation 
demonstrates robustness of the L; so sample adopted above. 




Figure 1. The £i so distribution and the fluence limit. 
3 STATISTICAL TECHNIQUES 

The E is0 - z sample suffers from various selection effects 
(Nava etal ■l2008l ; lLu et al.l201 1), among which the dominated one 
is the truncation due to the detection limit of the telescope, as seen 
in Fig. 1. If this bias not removed, the E iso — z correlation would be 
far from the intrinsic one. A nonparametric r statistical technique 
m ay be introduced to r esolve this problem, which was first put forth 
bvlLvnden-Belll (1 197 ll) and further d eveloped bv lEfron & Petrosiar] 
( 1992). Llovd-Ronnin g et all d2002h fir st applied this techni que to 
GRBs with simulated redshifts and later lYonetoku et all d2004b ap- 
plied it to a larger GRB sample still with simulated redshifts but up 
to z ~ 10. In this work, we continue to use this technique, but for the 
first time to two GRB samples with observed/measured redshifts. In 
short, this nonparametric t statistical technique uses a well-defined 
truncation criterion to estimate the correlation (if any) between the 
relevant variables and their underlying parent distributions. 

The £j so and z are not independent. Without loss of general- 
ity, the total isotropic-equivalent-energy function can be rewritten 
as ®(E lso ,z) = p(z)(f>(E lso /g(z))/g(z), where p(z) is the GRB forma- 
tion rate at z, <f>(E, m / g(z)) is the present-day (i.e., z = 0) isotropic- 
equivalent-energy function, and g(z) counts for the cosmic evolu- 
tion of E isa , that said, E'. = E iso /g(z). In the following analysis, the 
evolution g(z) will be removed from the E lso sample; after that one 
obtains the E' distribution, then the cumulative function il/(E' ), 

ISO ' v ISO ' 

and finally the GRB formation rate p(z). 

Consider a set of observable Ei so j and z, , where i indexes the 
zth burst and in our case i runs from 1 to 95. As shown in Fig. 1, for 
the ith sample of (z ; , -Eiso,;), we consider an associated set of 

Ji = U\E iso j > E iso ,i,Zj < z,,iiml, for 1 < < 95. (2) 

in which the number of samples in the 7, set is iV,-. The Zi,hm is the 
redshift of the crossing point between two lines of E = E iso j and 
the fluence limit corresponding to its "isotropic-equivalent-energy" 
limit. If z, and Ei so j are independent to each other, one would expect 
the number of the following sample 

Ri = Number^ e Mzj < Z;} (3) 

to be uniformly distributed between 1 and Np To estimate the corre- 
lation degree between E lso and z, one may introduce the test statistic 
t parameterized as 

Z,(Ri - Ed 



(4) 
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Figure 2. The statistical r as a function of the evolution parameter k. 
The best-fit r = as well as its lcr error (i.e., r = ±1) correspond to 
k = 1.80+jJ- 3 *. Thus, g k (z) = (1 + z) 1 80 is the best function to describe 
the isotropic energy evolution. A hypothesis of no evolution (equivalent to 
k = 0) is rejected at a significance of 2.5o". 



where = (Nj + l)/2 and V, = (A', 2 - 1)/12 are the expected mean 
and the variance of the uniform distribution, respectively. 

If Rj follows an ideal uniform distribution, then the samples of 
Rj < Ej and Rj > £, should be equal, and thus the statistic parameter 
t tends to be zero. Note that the r value here has been normalized 
by the square root of variance, so the correlation degree z and E iso 
can be measured in units of standard deviation. 

The two variables £ iso and z are connected through the 
form of E iso = 4ixT D 2 L (z) I '(1 + z), where T is the fluence in 
1-10000 keV in the burst restframe. Note that bursts in the sam- 
ple were actually detected by instruments onboard different satel- 
lites. The flux sensitivities of these instruments are as follows: 
~4x lO^ergcrrrV for CGROfBKTSE; ~ 10~ 7 ergcnr 2 s~ ] for 
Konus-Wind, BeppoSAX and Fermi/GBM; ~ 3 x lO^ergcirr 2 s -1 



for HETE-2; and ~ 10 erg cm 



for SwiftfBAT. Some bursts 



were detected by more than one satellite. To be safe, we choose the 
best sensitivity of ~ 10~ 8 ergctrr 2 s _I as the flux limit of the whole 
sample. 

For the Lj so sample, the two variables are connected by L- KO = 
4nFD^(z) where F is the flux in 1-10000 keV in the burst rest- 
frame. We set the flux limit as F = 5 x 10~ 8 ergctrr 2 s -1 in view 
of the fact that the sen sitivity of BAT is a few ~ 10~ 8 ergcirT 2 s _I 
( iBarthelmv et al.ll2005h . Nevertheless, our results are not sensitive 
to the limit of T and/or F. 



4 RESULTS 



Follow ing iMalonev & Petrosianl dl999h and ILlovd-Ronning et"al] 
(2002), we take the form of gt(z) = (1 +z) k in order to separate the 
isotropic energy evolution gk(z) from the GRB sample. The value 
E' iso = E[ s0 /gk(z) represents the isotropic energy after removing the 
evolution effect. When t is not equal to zero, we change the k values 
until t = with a proper k. Fig. 2 shows the r value as a function 
of k. The null hypothesis of the evolution is rejected at about 3.5 cr 
confidence level. The best fit to the E lso - z data yields that E' iso = 
E iso /a + z) ,S0 ,i.e.,k=l.S0. 



After converting £i s0 into E'. 



£ iso /(l + zV 



Figure 3. Cumulative isotropic-equivalent-energy function i//(E' iso ) of 
= £"150/(1 + z) 1 80 , which is normalized to unity at the dimmest point. 



equ ivalent-energy func t ion" iA(£,',„,) with the fo l lowing equation 
of jLvnden-Belll I1971I: lEfron & Petrosianl 1 19921 : IPetrosianl [l993l: 
ILlovd-Ronning et alj2002l : lYonetoku et alj2004l) 



ln^ 1 ' so ,,) = 2 ln ( 1 + ^-) 



(5) 



As can bee seen, the cumulative number at the rth point is 
calculated from Nj and for each point indexed by j, a truncation 
parallel to the axes is made and a weight l/Nj, based on the number 
of points in the associated set, is assigned to that data point. 

Fig. 3 shows the "isotropic-equivalent-energy function" of 
E' lso = £j so /(l+z) L8 °. The shape of the "isotropic-equivalent-energy 
function" roughly follows a broken power law, and the dim and 
bright segments can be parameterized as if/(E\ so ) oc Er° 21±om and 
'f'(Eiso) K ^o' 87 * 07 • It corresponds to the isotropic energy distri- 
bution at z = since the evolution effect has been removed. The 
"isotropic-equivalent-energy function" in the comoving frame is 
t«£' so )(l+ Z ) 180 . 

To estimate the cosmic GRB formation rate, the cumulative 
number distribution tf/(z) as a function of z is derived using the func- 
tion analogous to equation (5). In this case, for the rth sample the 
associated set is given by 



— ij\Zj ^ Zh E[s j > ^iso,(,lim}) 



(6) 



where E iso j ] im is calculated at the crossing point of the fluence limit 
and z = Zi- The resulting cumulative GRB formation rate if/(z) is 
shown in Fig. 4. 

To be scientifically useful, one needs to convert the cumulative 
formation rate into the differential form (e.g., to compare with the 
cosmic star formation rate). The conversion is given by 



dz 



dV(z) j 
dz J ' 



(7) 



where (1 + z) comes from the cosmological time dilation, and 
dV(z)/dz is a differential comoving volume described by 



dV 

lz 



H o Jo Jn, + n,„z 3 



1 



(8) 



can nonparametrically derive the cumulative (local) "isotropic- 



Fig. 5 shows the resulting differential GRB formation rate 
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Figure 4. Cumulative GRB formation rate tp(z) as a function of z, which is 
also normalized to unity at the highest point. 



Figure 6. Cumulative luminosity function ip{L[ so ) of V- = Li so /(l + z) 2,30 , 
which is normalized to unity at the dimmest point. The evolution effect is 
removed. 





Figure 5. Relative GRB formation rate, which is normalized to unity at the 
first point. The error bar represents the 1 <x statistical uncertainty of each 
point. 



Figure 7. Cumulative GRB formation rate of Lj so sample as a function of z, 
which is normalized to unity at the highest point. 



The best-fit power-laws for different segments are 

( (1 +z ) 8 - 24 ± 448 forz< 1 
p(z) °c \ (1 + z )-°- 54± ° 64 for 1 < z < 3.5 (9) 
( (i +z )-3.80*zi6 for z > 3.5 

with 95% confidence bounds. 

The L iso - z sample was treated in the same way. For this sam- 
ple, we found k = 2.30*°^, which is close to the value from the 
£j so - z sample. The corresponding luminosity function, cumula- 
tive GRB formation rate and differential form of the GRB forma- 
tion rate have been reported in Figs. 6-8. 

We use the best-fit isotropic-equivalent-energy function and 
GRB formation rate to calculate log N-log S distributions and com- 
pare them with the observed log N-log S distributions of BATSE 
and Swift bursts. The results are shown in Fig. 9. From this figure, 
we find that the results obtained from our model can well repro- 
duce the observation (The probabilities of the K-S test are P=0.95 
and P=0.74 for the BATSE and Swift GRB samples, respectively). 




I I I I I I L 

1 2 3 4 5 6 7 8 

1+z 



Figure 8. Relative GRB formation rate of L lso sample, which is normalized 
to unity at the first point. The err bar represents the 1 tr statistical uncertainty 
of each point. 
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Hopkins & Beacom (2006) 
-GRB rate: This work 

LGB: Mannucci et al. (2007) 
Verma et al. (2007) 
Bouwens et al. (2008) 

LAE: Ota et al. (2008) 



log(S) 




log(S) 



Figure 9. Cumulative bursts distribution as a function of the fluence (S). 
Panel (a): log N-log S distribution for BATSE bursts (Solid line) and the 
distribution predicted by our model (Dotted line). Panel (b): log N-log S 
distribution for Swift bursts (Solid line) and the distribution predicted by 
our model (Dotted line). 



5 CONCLUSION AND DISCUSSION 

GRBs are brief but intense emissions of soft y-ray, lasting from a 
few seconds to a few thousand seconds. For GRB-like high energy 
transients, the total isotropic-equivalent energy, £ iso , can be reli- 
ably measured, and the number density of bursts per £; so interval 
may provide an independent or even more representative clue on 
the underlying physics of GRBs. In this work, using a sample con- 
taining 95 bursts with measured redshifts, we for the first time con- 
structed the isotropic-equivalent-energy function and then adopted 
it to estimate the GRB formation rate. The fluence-truncation ef- 
fect has been properly addressed by adopting a t statistical tech- 
nique. We find there exists cosmic evolution between £i so and z, 



i.e., E im oc g k ( z ) = (1 + z) 



(see Fig. 2), which is comparable 



with that between L iso and z, i.e., L iso oc (1 + z) . Our finding 
of the evolution is largely consistent with previous findings using 
GRB samples with simulated redshifts jLlovd-Ronning et al . 2002; 
IWei & Gadl2003l ; lYonetoku et al.ll2004t) . 

We also find that the evolution power index k is not sensitive 
to the chosen value of fluence/flux as long as it is around the instru- 
mental sensitivity. Changing the fluence limit from 8x 10~ 7 erg cirr 2 
to 5 x 10~ 8 ergcirT 2 , the best-fit index k of E ls0 sample only varies 
from 1.80 to 2.05, changing the flux limit from 5 x 10~ 8 erg s~'cm~ 2 



Figure 10. GRB ra t e com p aring with SFR. Th e data of SFR come from 
iHopkins & Beacorrt fe00d) ; lKistler et all l2008l) ; lYiiksel et all f200l . We 
normalized the first point of GRB rate and SFR for convenience. 



to 1 x 10~ 8 erg s~'cirr 2 , the best-fit index k of L iso sample only varies 
from 2.30 to 2.65, all well within the lcr uncertainty. 

As can been seen, the evolution of £i so is comparable with 
that of L iso . This might indicate that the durations of GRBs do 
not evolve significantly with redshift and that the existence of the 
isotropic-equivalent-energy evolution may reflect the evolution of 
typical physical parameters of GRB progenitors. 

After removing the redshift dependence, the isotropic- 
equivalent-energy function can be reasonably fitted by a broken 
power-law. For the dim and bright segments, we have if/(E isa ) oc 



and if/(E ISO ) oc E. 



respectively (see Fig. 3). The 



shape of the energy function (see Fig. 3) is similar to that of 
the luminosity function (see Fig. 6). Moreover, our indices of 
(-0.27, - 0.87) are comparable t o the indices of (-0.29, - 1.02) 
reported bv lYonetoku et alj d2004l) . 

The GRB formation rate as a function of redshift is 
also calculated. The connection between (lon g-duration) GRBs 
with broad-lined Type Ic supernovae (e.g . . IStanek et alj 20031 ; 
iHiorth et alj|2003l ; IWooslev &BloorrJ 120061 ; iFan et alj|201lh sug- 
gests that GRB progenitors are very massive, short-lived stars, 
leading to expectation that the cosmic GRB formati on rate would 
nicely follow the c osmic star format i on hi st ory l Totanil [79971; 
I wiiers et all Il998l; ILamb & Reichartl l200d ; iBlain & Nataraiarj 
20001; IPorciani & Madaul200lh . So far, the observed star formation 
rate follows p(z) oc (1 + z) 3 4 for z < 1, oc (1 + zT 3 for 1 < z < 4 
and oc (1 + z )- 3 5 for z > 4 dKistler et al]|2008l;lYuksel et al.ll2008h 
(see Fig. 10). Our results (see Fig. 10) show that the GRB forma- 
tion rate increases quickly for z ^ 1, then roughly keeps constant 
for 1 < z < 4, and finally decreases at higher redshift with a power 
index of -3.8 (i.e., oc (1 + z)~ 3 ' 8±2 ' 16 ), in good agreement with the 
star formation rate. In future a larger GRB sample with measured 
redshifts would better address this topic and make progress in the 
fields of massive star formation and GRB physics in the early uni- 
verse. 

Our results are based on the assumption of a power-law red- 
shift evolution of the £i so . Besides introducing the evolution of the 
isotropic-equivalent-energy function, other models can also repro- 
duce the observed data, such as one could enhance high-z GRB rate 
by introducing a metallicity preference of GRBs or increase the to- 
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tal number of GRBs (e.g., Virgili et al. 201 1). These models might 
be distinguished by the future observations. 
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